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We compute strain distributions in core-shell nanowires of zinc blende structure. We use both 
continuum elasticity theory and an atomistic model, and consider both finite and infinite wires. The 
atomistic valence force-field (VFF) model has only few assumptions. But it is less computationally 
efficient than the finite- element (FEM) continuum elasticity model. The generic properties of the 
strain distributions in core-shell nanowires obtained based on the two models agree well. This 
agreement indicates that although the calculations based on the VFF model are computationally 
feasible in many cases, the continuum elasticity theory suffices to describe the strain distributions 
in large core-shell nanowire structures. We find that the obtained strain distributions for infinite 
wires are excellent approximations to the strain distributions in finite wires, except in the regions 
close to the ends. Thus, our most computationally efficient model, the finite-element continuum 
elasticity model developed for infinite wires, is sufficient, unless edge effects are important. We 
give a comprehensive discussion of strain profiles. We find that the hydrostatic strain in the core 
is dominated by the axial strain-component, Ezz- We also find that although the individual strain 
components have a complex structure, the hydrostatic strain shows a much simpler structure. All in- 
plane strain components are of similar magnitude. The non-planar off-diagonal strain-components 
(sxz and eyz) are small but nonvanishing. Thus the material is not only stretched and compressed 
but also warped. The models used can be extended for study of wurtzite nanowire structures, as 
well as nanowires with multiple shells. 

PACS numbers: 62.23.Hj, 62.35.-g, 62.20.D-, 68.70.+W 



I. INTRODUCTION 

Nanowires represent a promising technological platform with a wide range of applications spanning from electronic^ 
and photonic deviced^HH to biochemical sensor^. The reasons for this success are not only the increasing ability in 
handling and manipulating nanowires but also an ongoing improvement of the quality of the crystal structure of the 
grown wires. 

The progress in crystal growth has led to the fabrication of coherent crystalline nanowires with core-shell structure^. 
The differences between the materials of the core and shell in lattice constants and band parameters allow for tuning 
the properties of the resulting wire^. Possible strategies involve confining charge carriers to the nanowire core to 
reduce the influence of the surface on the electrical properties. Likewise, a core-shell nanowire can serve as an 
optical waveguide confining light. Confining charge carriers to the shell, or having multiple active shells and confining 
electrons and holes to different shellaSEl, may also be relevant, depending on applications. The confinement can be 
obtained by band gap engineering and strain-engineering^. In strain-engineering, one chooses to grow the core and 
shell in a nanowire with different materials with a coherent crystalline core-shell interface. Thus the nanowire is under 
a pseudomorphic strain, due to the deformations at the interface to accommodate the different lattice constants. This 
strain generates different deformation potentials in the core and shell regions which affect the electronic structure 
and thus the electrical and optical properties of the nanowires. In particular, gaining deep insight into the electronic 
structure of heterostructured nanowires with a lattice mismatch, requires finding the underlying elastic deformation, 
since this elastic information is needed as input to, e.g., k ■ p or tight-binding calculations.^"^ 

Theoretical investigations based on the continuum elasticity of systems experiencing pseudomorphic strain go back 
to Eshelby^ who used analytical methods to study the strain effects of one material immersed in another. More 
recently, several studies of pseudomorphic strain fields have been performed for finite nanostructuresP^ On core-shell 
nanowire geom etries, Niquet studied the effect of a shell on the electronic properties of a quantum dot embedded 
in a wireP2ED Studies on the effect of a core-shell geometry on the electronic properties in quantum wires were also 
performed by Pistol and Pryorpl Schrier et a^and Zhang et al. 9 However, these studies only considered very small 
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structures or presented the hydrostatic strain in the structures. 

In this article, we calculate and discuss the elastic deformation in core-shell nanowires of zinc blende structure 
with a lattice mismatch between the core and shell materials. Numerical calculations are performed for both finite 
and infinite wires using both continuum elasticity theory (implemented by a finite clement scheme) and the atomistic 
model of Keating 2 -^. We present the strain distributions in a cross section of the nanowire as well as line scans 
both in the cross-section and along the wire. We consider free nanowires, i.e., we neglect external forces acting 
on the nanowires. Core-shell nanowires with various cross sections are investigated: hexagonal cross sections with 
both parallel and non-parallel core and shell facets and co-centric circular cross sections. We restrict ourselves to 
only discussing nanowires with a single shell and a core, although the methodology used here can be generalized to 
multiple shell nanowires. Likewise, we shall only consider purely zinc blende crystalline heterostructures. 

The article is organized as follows. The finite-element continuum elasticity model employed in this work is described 
in Sec. [II] and the atomistic model in Sec. |III| Numerical results and discussion are presented in Sec. |IV| Finally, Sec.[V| 
is devoted to summary and conclusions. 

II. CONTINUUM ELASTICITY THEORY 

We here describe the continuum elasticity (CE) model employed for a core-shell nanowire of zinc blende structure 
with axis along the [111] direction for which the material lattice parameters vary discontinuously at the core-shell 
interface. 



A. Coordinate systems 




FIG. 1: Nanowire with core and shell. Our results refer to a core of GaAs inside a shell of GaP. The core hexagon has a 
side-length of 6 nm and the shell hexagon has a side-length of 13.9 nm. The axis direction is [111]. 

We consider core-shell nanowires of zinc blende structure with axis along the [11 Indirection. Figure [I] shows a 
typical core-shell nanowire geometry. We define two coordinate systems: a standard coordinate system with axes 
denoted by (x,y,z) coinciding with the main crystallographic axes ([100], [010], [001]) and another system where the 
axes will be denoted by (X, Y, Z), referred to as the nanowire-oriented coordinate system. In the (X,Y,Z) system, 
the Z axis is parallel to the axis direction of the nanowire (i.e., the [111] direction). The X and Y axes are in the [110] 
and [112] crystal directions, respectively, as depicted in Fig. ^1 Formally, the nanowire oriented coordinate system is 
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FIG. 2: Cross-section view of the core-shell geometry. Coordinate axes X and Y indicate our nanowire oriented coordinate 
system. The indicated paths A and B are used when we present our numerical results in Sec. |IV| 



related to the standard coordinate system by the transformation 
with a rotation matrix 1Z given by 

K = 



(1) 




(2) 



The strain energy U is written as 



B. Strain energy 



U = JwdV=^ J Cjjki Ejj gfcj dV 

ijkl 



(3) 



where w is the strain energy density, Cijki are the elastic stiffness tensor elements, and Sij are strain tensor 
componentiP^. 

For a material with cubic symmetry the elastic stiffness tensor, cyfcj, only has three independent elements and the 
strain energy density in the crystallographic coordinate system can be written as 

w = \ c n{ £ lx + £ lv + £2 zz) + ci2{e xx e yy + e xx e zz + e yv e zz ) + 2c 44 (e 2 xy + e 2 xz + e 2 yz ), (4) 



where 



Cl2 = C%. 



C44 



(i = x,y,z) , 

"ijji {hj = x, y, z, i ^ j) 



(5) 



The identities in Eq. ^ reflect the cubic symmetry of the material and define the general symmetry properties of the 
elasticity tensor.^ In our calculations, numerical values of the are taken from Ref. [251 



C. Incorporating pseudomorphic strain 



In this section we describe the lattice matching conditions and their interpretation in two- and three-dimensional 
bulk systems. 
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1. Pseudomorphic conditions 

We consider two different elastic media, which are the core and the shell materials of the same crystal structure. 
We shall use superscripts (c) and (s) to indicate core and shell, respectively. A typical geometry for a finite wire is 
shown in Fig. [T] and a cross section of the wire is shown in Fig. [2] 

For cubic materials the lattice vectors a( c ) and a( s ) of core and shell in the strain-free state can be written as 



a 



to =\^to 



5>toe, = ato£<, 



aW=^a W e i =ato^n i e i , (6) 



where and ato are the lattice constants of the core and shell materials, af' and a> are the components of the 
lattice vectors in the core and shell, n., are integers, and are the basic vector in the crystallographic coordinate 
system. 

After lattice deformation, distorted lattice vectors in the core and shell can be written in terms of the components 
efj and e\j of the corresponding strain tensor, as 



y J 


(c) 




)a) 



(7) 

At the core-shell interface, the in-plane component of the distorted vector r( c ) on the core side must match its 
equivalent r^ s ^ on the shell side. This pseudomorphic requirement can be expressed as 

E^[^ +£ S ) ] a 5 c)= E^h +£ «] a i s) ' (8) 

ij ij 

where ij are the components of an arbitrary tangent vector t of the interface. One implication of Eq. p| is that the 
strain will be discontinuous at an interface between two cubic materials having different lattice constants. 



2. Initial strain 



To fulfill the conditions of Eq. Q, we assume the strain tensor in Eq. ^ is given by a sum of two terms 



e= - [V ® u + (V ® u)*] + e (0) , (9) 



where u is the displacement field relative to a matched (yet unrelaxed) configuration and e^ ^ is an initial strain^. 
The initial strain is assumed nonvanishing only in the shell, where it is defined, in the cubic case, by 

^ = "—^.1 (10) 

where and are the lengths of the lattice vectors in the core and the shell, respectively. This initial strain 
ensures that the conditions Eq. ([8| are satisfied and leads to a matched structure. 

The structure obtained for u = is not that of lowest energy for =^ and the structure is relaxed by varying 
the field u to minimize the energy Eq. (|3]) . For a finite wire, the degrees of freedom to vary are precisely the degrees 
of freedom of u. For the case of an infinite wire, the periodicity in the Z-direction is utilized and the resulting model 
is described below. 



D. Infinite Wire 



The model of an infinite wire utilizes the periodicity in the Z-direction and reduces the modeling domain to a 
2-dimensional domain. 
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The strain tensor is transformed to the nanowire oriented coordinate system by 



£ e>!/3 = ^ao^li e ab ■ (H) 
ah 

where Greek indices vary over X, Y, and Z and Latin indices vary over x, y, and z. 

In the model of an infinite wire, the vectorial field u introduced in Eq. ^ depends only on the coordinates in an 
(X, Y) plane of the nanowire, and we write u(X,Y) G M 3 . The axial strain in the core e zz an d the axial strain in 
the shell e zz are constant and by Eq. (JsJ) fulfill 



1 + efz 



with the constants 



a^ s > . (12) 
The axial strain can thus be described as a single degree of freedom a as 

= — — 1 with i = s or c . (13) 

zz am v ' 

In the model of an infinite wire, the field u(X, Y) and the variable a are the degrees of freedom. 
The energy density [defined in Eq. ([3])] becomes 

w = 7,[ F i( £ xx + £ yy) + f 2£xx£yy + F 3 (e X x + £yy)szz 
+F±{exx — £yy)zyz + ^F^exz^xy + F§e zz 

+F 6 (e YZ +e 2 xz ) + F 7 e XY ] (14) 

Fx = * (cu + cia + 2c 44 ) , (15) 

^2 = ^ (en + 5ci2 - 2c 44 ) , 
2 

^3=2 (cu + 2ci2 - 2c 44 ) , 
-F4 = — ^— (en - C12 - 2c 44 ) , 

^5 = 3 (cu + 2ci 2 + 4c 44 ) , 
4 

-F6 = g (en - C12 + c 44 ) , 
2 

Fr = - (en - C12 + 4c 44 ) . 

In this model, w e keep £xz and £yz in contrast to the conventional plane strain approximation in which exz = 
£yz = is assumed 22122 '. 

The minimization of the strain energy over the variable a corresponds to the condition that the total axial force Fz 
vanishes. This has been discussed for a simpler energy functional in Ref. The vanishing of the total axial force 
can be understood by considering a wire of a fixed number of unit cells N. By translation invariance the total energy 
is proportional to N. The variation of the energy per unit cell becomes 



1_ (dU c dU s \ 
N\da da) 



where U c and U s are the strain energies of core and shell. The energy is given by Eq. ([3| with the integration taken 
over the undeformed domains 2 ^. For the core the undeformed axial length l c is proportional to Na^ and for the shell 
the axial length l s is proportional to Na^ s ' . In particular, the total energy of the core with energy density w c is given 
by 

U c = [ dz ( dSw c <x 7Va (c) [ dS w c . (17) 
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For the variation of the core energy, we calculate the partial derivative 

dw c _ 1 dw c _ a Z z 
da a( c ) dezz 



(18) 



where azz = F^ezz + F^(exx +£yy)/2 is the axial stress. The area element dS and the axial stress then correspond 
to an infinitesimal axial force dF z = dS azz- Using these results and those similar for the shell, we can write 



1 dU 

N ~da 

oc fdSazz + I dSazz = I dF z + I dF z = F z 



"=NlTa (19) 



proving the assertion. This result [Eq. ( 19 1] could be generalized to infinite wires under stretch. In that case we would 
use a non-zero force condition. 

The actual minimization is done using finite elements^ (linear triangular elements) for the continuous displacement 



u and a variable a to parameterize the axial strain [Eq. (13)]. The finite element method (FEM) is chosen due to its 



flexibility with respect to arbitrary geometries. The method leads to a sparse matrix problem of the form F = K x 
from which x (representing u and a) is solved by standard numerical packages. Note that here the variable a couples 
to all other degrees of freedom in the matrix K. This is unlike the usual setting of only local coupling in two- and 
three dimensional bulk structures. Still, in practice the solution is very fast even for meshes with around 30000 nodes 
and three degrees of freedom per node. 



E. Finite Wire 



We also compute the strain in finite core-shell nanowires. The strain is computed by minimizing the total strain 
energy, defined in Eq. [3] This is performed in the coordinate system x, y, and z of the principal axes of the underlying 
crystal and by incorporating the pseudomorphic strain as described in Sec. II C For finite wires we have no benefit 
of using the nanowire oriented coordinate system because of the lack of translational symmetry. 

The typical length L z of a modeled finite wire range from 30 nm to 150 nm. In the modeling of finite wires we 
exploit the C^ v symmetry of the nanowire. That is, we model and mesh only a segment of 1/6 of the total nanowire 
geometry. On the nodes on the side facets of this segment we then fix the displacement perpendicular to the facetP^. 
By this symmetry reduction we reduce the computational domain and the number of nodes by a factor of 6. This 
technique also reduce the source of numerical errors since the resulting element mesh has the same symmetry as the 
nanowire geometry. 

The CE model of finite wires is implemented using three-dimensional tetrahedron elements with second order 
polynomials. A typical element mesh contain in total around 1.5 • 10 5 nodes distributed in a wire, with 3 degrees of 
freedom per node. We use a nonuniform and adapted element mesh with at most 9 • 10 4 nodes per cross-sectional 
area (on 1 /6 of the nanowire cross-section) . 



III. ATOMISTIC VFF MODEL 

For the atomistic calculations we use the Valence-Force Field (VFF) model of Keating^!. The structure is built up 
one atomic layer at a time. We choose a layer structure that corresponds to a zinc blende phase nanowire with axis 
along the [111] direction. 

We first describe our model for finite wires and then describe our implementation for an infinite wire using translation 
invariance. 



A. Finite wires 



The atomic coordinates {ri}fL 1 are the degrees of freedom of our VFF model. The energy is a function of these 
and is written as sums over interatomic bond distances and three-body bond angles. 
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We let Tij = Vj — Ti denote the vector from atom i to atom j and r^o denote the value ry takes in a stress-free 
bulk crystal. The potential energy of the crystal is given by^ 

E = EE^(4-4o) 2 (20) 



where the equilibrium bond angles fulfill r ij0 ' r iko — ~ \ r ijo\ l r ifeol /3 in zinc blende and n.n. on top of a summation 
sign indicates that we sum over nearest neighbors of atom i only. 

Numerical values for the coupling constants a and f3 are obtained by comparing the energy in the VFF model 
Eq. ([20} and the CE model Eq. and requiring that the energies coincide for small deformations of a small volume. 
This requires three equations to be satisfied, one for each material parameter in the CE model. As the VFF model only 
contains the two coupling constants a and /3, all three equations cannot be exactly satisfied, but any choice of values for 
a and j3 will give rise to effective elastic constants of the VFF model that deviate form the experimentally determined 
values used in the CE model. We follow Pryor et. alP^ in using the values provided by Martin 31 . Numerical values 
for the coupling constants and the relative errors in effective elastic constants are summarized in Tab. [Ij 



TABLE I: Coupling constants used in the VFF model, and relative errors for effective elastic constants. 



GaAs GaP 

a (N/m) 41.19 47.32 

P (N/m) 8.95 10.44 

r (A) 2.448 2.36 

Sen (%) -0.22 2.7 

5c 12 (%) 1.7 9.1 

<5c 44 (%) -13 -11 



Possible extensio ns of our model include an additional Coulomb interaction energy^ and higher order terms for 
non-linear elasticit y^ 2 ! 33 !. Such extensions are beyond the scope of the current work. The current model is widely 
applicable for analyzing linear elasticity of, e.g., zinc blende and wurtzite crystals. It could also be applied to the 
description of epitaxially strained heterointerfaces between materials with different lattice structures (see e.g. [3]). It 
is, however, worth to note that the applicability of the VFF model is, in this case, dependent on a priori information 
about the interface structure. 



B. Infinite wires 



For the simulation of an infinite wire, we start with a finite wire segment. We consider this finite segment to be a 
unit cell of an infinite wire consisting of an infinite array of such unit cells. The distance between two unit cells in the 
axial direction is added to the model as an additional degree of freedom Ar and is called the axial lattice constant in 
the following. 

Bonds between the different unit cells are added to the energy expression for a finite wire segment which was given 
in Eq. ( 20 ) . This gives an expression for the energy per unit cell of an infinite wire. The inter-cell bonds couple the 
degrees of freedom of the topmost and bottommost atoms, via the axial lattice constant Ar. We then minimize the 
energy per unit cell, the degrees of freedom now being the atomic positions inside the unit cell and the axial lattice 
constant Ar. 



C. Definition of strain 



To enable comparisons with continuum elasticity we must define a deformation measure from the atomistic VFF 
displacements that is comparable to the macroscopic elastic strain. Strain is a macroscopic concept and several 
different atomistic strain quantities can be defined, that all converge to the macroscopic strain in the relevant limit. 
In this work, we follow the strain definition of Pryor and coworkers^. 



t 



FIG. 3: The tetrahedron of neighbors around one atom in the atomistic model. 



Microscopic definition 



The microscopic definition of strain at the position of an atom is formulated in terms of displacements of the atoms 
in its local neighborhood. For an atom, positioned at ro, the relative position vectors r%, . . . of its neighbors form 
a tetrahedron around it as shown in Fig. [3} In order to define our concept of strain at r , we consider these 5 atoms 
to form the local neighborhood of the atom at ro. Our definition of strain will use a deformed and an undeformed 
configuration of this neighborhood. In the deformed configuration, the positions of the 5 atoms are the same as in 
the result of the full VFF calculation. In the other, the undeformed configuration, the positions of the atoms are 
chosen such that the energy of the isolated 5-atom system considered, as defined by the VFF model (20), vanishes. 



In the cases of a pure GaAs or a pure GaP tetrahedron, the undeformed configuration is equal to a stress-free bulk 
configuration. At interfaces, there is no corresponding bulk system. 

For each of the two configurations, we define the vectors r.y = rj — r; from the coordinates of the tetrahedron 
corners, and from these vectors we construct a 3 x 3 matrix whose columns are given by r 2 i, r 32 and r^. Strain is 
then definecPSl via 



_(<*) (d) (d) 

l 21 l 32 l 43 



(1 

£ - 



,(0) JO) (0) 



L 21 



L 32 



1 43 



(21) 



where the superscripts d and 0, refer to the coordinates of the deformed and undeformed state, respectively. The 
strain at ro is, consequently, the symmetric part of the tensor that describes, to lowest order, the deformation of a 
small volume around the atom at Tq. 



2. Numerical average 



For the microscopic strain defined by (21), we find strong oscillations in the strain field. Neighboring atoms can 
have very different strains. This is particularly clear at the GaAs-GaP interface, where the microscopic description of 
the interface renders it non-planar (even locally), and this is phenomenologically different from the continuum elastic 
description of the same interface. In order to allow for a comparison of the strain, obtained with these two elasticity 
models, we have chosen to smoothen the strain of the VFF model. We perform this smoothing by averaging the strain 
over each pair of group III and group V atoms aligned along the [111] direction. 

The static strain field of continuum elasticity describes macroscopic features of the continuous material, and thus 
implicitly contains a smoothing over any internal deformation at the atomic scale. The strong variations we find on 
this scale can thus have no direct correspondence in continuum elasticity, and our averaging procedure removes those 
effects. 



IV. RESULTS 



When comparing the results for finite and infinite nanowires, we evaluate the strain of a finite wire at cross sections 
in the VF-plane far enough from the free ends of the nanowire. We only discuss the strains in the (X, Y, Z)-system. By 
the principle of Saint- VenanlP^l, the actual force and moment distributions at the cross-sections are not important for 
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X(nm) X(nm) X (nm) 



FIG. 4: Strain in an axial unit cell of an infinite nanowire with a hexagonal geometry. The results were obtained with the infinite 
atomistic model. The plots (a) to (g) show the strain components exx, £yy, £zz, £yz, £xz, and £xy and the hydrostatic 
strain £xx + £yy + £zz, respectively. The axes are as indicated in Fig. [2] 
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FIG. 5: Strain fields in a unit cell of an infinite nanowire where the core and shell hexagons have the same orientations. Results 
obtained from the infinite atomistic model. The plots (a)-(g) show the strain components £xx, £yy, £zz, £yz, £xz, and £xy 
and the hydrostatic strain £xx + £yy + £zz, respectively. 

the state far enough from the cross sections. Only the total force and moment matter. For a meaningful comparison, 
we therefore require that the resultant force and moment on these ends are similar to the conditions on the cross 
section of the infinite wire. In our case, we consider free nanowires without external forces. This corresponds to no 
total force on the cross-sections. As we expect, the results for a cross section of the finite models converge to those of 
the infinite models with increasing distance of the cross section from the free ends. 

Deviations between the atomistic and continuum models are expected to be dominated by the difference in effective 
elastic constants. The force parameterization in the VFF model uses only two parameters that are to be determined 
from the three material parameters of cubically symmetric continuum elasticity. We find that this imperfect matching 
is indeed the main source of discrepancies between the results of the VFF and CE calculations. 

A. Strain distribution 

Figures [4] to [6] show the strain fields in the X V-plane for an infinite wire modeled using the VFF model and three 
different core-shell geometries. The strain of the hexagonal geometry (depicted in Figs, fl] and [2]) is displayed in Fig. El 
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FIG. 6: Strain in a unit cell of a nanowire with cylindrical core-shell geometry. Results obtained from the infinite atomistic 
model. The plots (a)-(g) show the strain components Exx, £yy, £zz, £yz, £xz, and £xy and the hydrostatic strain £xx + 
£yy + £zz, respectively. 



All strains are given with respect to the nanowire oriented coordinate system (X, Y, Z) . Corresponding plots for a 
parallel hexagonal geometry, where the core hexagon is rotated by 30° relative to that of Fig. [2] and for a cylindrical 
geometry, are shown in Figs. [5] and [6| In all cases the nanowire is composed of a GaAs core and a GaP shell. We 
discuss the strain components in the hexagonal geometries, and then compare with the cylindrical geometry. 



1. Axial Strain 



As discussed in Sec. |II C 1| we see in several of the figures that the strain is discontinuous at the core-shell interface. 
The axial strain ezz is displayed in Fig. |4]jc), and we see that it is constant within core and shell (as was used in 
the infinite CE model). The core material (GaAs) has a larger lattice constant than the shell material (GaP), and 
the energetically most favorable core-shell nanowire configuration will have an interatomic spacing that is in between 
those of bulk GaAs and bulk GaP. We therefore see an axial compression (ezz < 0) of the core and an axial tension 
(e Z z > 0) of the shell. 



2. Hydrostatic Strain 

We see from the hydrostatic strain in Fig. Bg) that the core is compressed while the shell is expanded. In the 
shell, the magnitude of the hydrostatic strain [Fig. [4^g)] is similar to the axial strain [Fig. [4^c)] . Thus in the shell, 
hydrostatic strain is dominated by the axial contribution. 



3. Planar Tensile Strains 



Tensile strains correspond to expansion or compression along the coordinate axes. In Fig. |4|a), exx, i-c, the 
compression in X-direction, is shown. The structure is compressed in the X direction everywhere along the X-axis. 
The core is compressed, as it has to fit into a space too small for it to be at equilibrium, and the shell is compressed 
along this axis because it is pushed outwards by the core. Along the y-axis, the core is again everywhere compressed, 
but the shell has been expanded relative to its equilibrium state. As discussed for the axial strain, the equilibrium 
configuration will, parallel to an interface, have an interatomic spacing in between those the two materials take in 
bulk. Therefore, the shell, with a smaller bulk lattice constant, is expanded. A trace of the same effect can be seen 
in the Y component of the strain, Fig. |4](b), near the two corners of the core with Y = 0. In this case, there is no 
interface exactly parallel to the y-axis. Nevertheless, the effect is still visible at the corners. In Fig.^b), the core 
has been rotated by 30° relative to the geometry of Fig. Qb). In this case, there is an edge of the core-shell interface 
parallel to the Y axis, and we see an expansion in the shell close to that edge. 
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4- In-plane Shear Strain 

The shear strains £xy, £xz and eyz are nonzero if the local deformation changes the angles between the basis 
vectors. Considering the lower left part of Fig.[4]jf), and applying the same argument we used for the strain components 
Sxx and eyy, we expect radial compression and circumferential extension of the shell, and we see a nonzero planar 
shear strain Exy- The negative value corresponds to the fact that the right angle between the unit vectors X and Y 
is locally enlarged by the deformation^. Analogously, in the upper left part of the same figure, the angle between X 
and Y is locally shrunk, and this corresponds to a positive value for £xy ■ 

5. Warp Effect 

The nonzero values for eyz and Sxz, shown in Fig. [4jd) and Fig. [4^e) , correspond to the fact that a cross section in 
the W-plane is warped into a non-flat surface in the strained wire. In other words, when the core pushes outwards 
on the shell, it is energetically favorable for the shell to respond by not only becoming compressed radially, but also 
to deform by warping out of the W-plane. 

6. Cylindrical Geometry 

Fig. [6] shows numerical results for a cylindrical nanowire geometry. The general features are similar to the previously 
described geometries. The symmetry of the strain due to the zinc blende lattice can be seen more clearly in the 
cylindrical geometry. Figures |6ja) and (b) show that exx and eyy are almost identical in the shell, except for the 
sign. In a significant part of the shell, they are larger in magnitude than both the axial [Fig. (6^c)] and the hydrostatic 
[Fig.J^f)] strains. We note also that the warp strains, shown in Fig.[6]Jd) and Fig.[6|e), do not vanish in the cylindrical 
wire. Thus, the warp effect is not solely an artefact of the hexagonal nanowire geometry, but an effect of the zinc 
blende lattice structure. 



B. Linescans 

In Fig. [7j strain components are plotted along paths which are parallel to the nanowire axis. The solid line 
corresponds to a path in the core [(X, Y) — (2nm, lnm)], and the dotted line corresponds to a path in the shell 
[(X,Y) = (8nm,4nm)]. 

As discussed previously, far away from the free ends, the strain converges towards its value in an infinite nanowire. 
We find that, for points that lie more than 30 nm away from the ends of our finite nanowire, no strain component 
deviates by more than 5 • 10~ 6 from its value at the center of the finite nanowire (Z = 0). We interpret this as an 
indication that our nanowire geometry has sufficient length and that the strain field at Z = will correspond to the 
strain field in the infinite model. 

The behavior of the strain as a function of Z can be understood as a combination of two effects. The first effect is 
due to the hexagonal core-shell geometry and the zinc blende lattice structure and was discussed above. The second, 
the Poisson effect, dominates the strain in the vicinity of the free ends. The core of the wire is radially compressed 
by the shell, and responds by expanding significantly along the axial direction. Far from the ends, this expansion is 
prevented by the shell, and the core is axially compressed, as seen in the plot of ezz in Fig. [7J Close to the free ends, 
the core bulges out of the wire, giving the end-surfaces of the nanowire slightly convex shapes. This allows the core 
to relax more, as is seen in the same plot, where the magnitude of ezz decreases towards the end of the wire. The 
Poisson effect is responsible for the main features of the Z-dependence of the other strain components as well. 

Figures [8] and [9] show a comparison between the strains of the different models along paths in an W-plane. The 
model not shown is the VFF model for a finite wire. Far from the ends, this model agrees well with the infinite VFF 
model. We observe that the results from the continuum models for finite and infinite wires are very similar. 

Deviations between the VFF and CE models at the core-shell interfaces are expected. These deviations are un- 
physical, in the sense that macroscopic strain is a locally defined phenomenon within a homogeneous material. We 
also see significant systematic deviations between the VFF and CE models inside the core, and we attribute these to 
the discrepancy in effective elastic constants of the VFF and CE models. As a test of this, we compared calculations 
from the CE and VFF models, where we tuned the parameters of the CE model so that the two models used the same 
effective clastic constants. A comparison of the strain fields from those calculations showed good agreement between 
the VFF and CE models (comparable to the agreement between our finite and infinite models). 



12 



0.0001 



-0.0001 



0.0001 

£ 

-0.0001 





N 

N. 

-0.0002 



0.0001 

N 

r 

-0.0001 
0.0002 



-0.0002 



i~ 

w 

-0.0002 

20 40 

Z (nm) 



FIG. 7: Linescans of the strain along the nanowire axis. The full-drawn lines correspond to (X,Y) = (2nm, lnm) (located in 
the core) and the dotted lines to (X, Y) = (8nm, 4nm) (in the shell). Results obtained from 3D continuum model. 

V. CONCLUSIONS 

We have given a self-contained comprehensive account of our computations of strain distributions in finite and 
infinite core-shell nanowires with lattice mismatch using both continuum-elasticity theory and an atomistic VFF 
model. The strain profiles were shown in Figs. [4] to [9] and were discussed. The atomistic VFF model has in this 
work been formulated within a unit cell of a nanowire with periodic boundary conditions in the axial direction. A 
corresponding continuum elasticity model for an infinite wire was also defined in a unit cell. Finally, a long and finite 
wire was studied using the continuum-elasticity theory and the VFF model. 

We observe bulging effects at the free ends of the finite wires and find good agreement between results at a cross 
section of a finite wire far away from its ends and at a cross-section of an infinite wire. The strain distributions obtained 
in the calculations using the continuum and atomistic models show qualitatively good agreement. Quantitatively, 
deviations remain, but their origin was discussed. In contrast to the hydrostatic strain, which has a very simple 
structure, the individual strain components have very complex structure. The non-planar shear strains are generally 
smaller than the diagonal components, but non- vanishing. The strain in the core is dominated by the axial strain, 
but at the core-shell interface and just outside it, a richer structure is found in the strain distributions. 

Future work will consist of computing properties such as electronic structure and phonons based on the results 
presented here. In those calculations, extensions to multi-shell structures and heterostructures combining zinc blende 
and wurtzite crystals will be considered. 
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FIG. 8: Linescan along path A as denned in Fig. [2] Results from three of the four models used are shown. The i-axis indicates 
position along the path. 
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